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It is demonstrated on a realistic model of amorphous alloy Sio.gGeo.i with 1000 atoms, that short- 
range spectral fluctuations of propagons and diffusons are universal and in agreement with ran- 
dom matrix theory. The universality ceases at distances greater than Thouless number Nt, where 
propagons obey the Altshuler-Shklovskii power law, £ 2 ~ N 3/2 for the variance £2 of the number 
of levels N, while a new power law, £2 ~ N, is observed for diffusons. 
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Complexity of the atomic structure of glasses results 
in a wonderful variety of vibrational eigenmodes [Q-Q- 
Topological defects and the lack of a long-range order 
lead to mode localization, inhibition of ballistic motion, 
or resonance; only a small portion of the modes resemble 
phonons-vibrations in periodic crystals. It is a difficult 
computational task to build a respectable-size atomic 
model of a glass and obtain the vibrational eigenstruc- 
ture by diagonalizing huge dynamical matrices; no an- 
alytical treatment is available. By applying the meth- 
ods of random matrix theory ||], this Letter reports on 
certain statistical analytical rules, universal (parameter 
free) and not, shared by similar modes. In particular, it 
is found that short-range spectral correlations of diffu- 
sons (non-propagating extended modes) and propagons 
(phonon-like modes) are universal and agree with random 
matrix theory, similarly to what have been observed in 
many other complex physical systems j^]. However, cor- 
relations at distances greater than Thouless number Nt 
cannot be described by random matrix theory. Spectral 
fluctuations become idiosyncratic: the statistics acquires 
a parameter (Nt) and the functional dependence is dif- 
ferent for different mode classes. Propagons follow the 
Altshuler-Shklovskii "3/2" power law for the fluctua- 
tions £2 of the number of levels, while diffusons have the 
fluctuations growing linearly with increasing number of 
levels. This is likely a general feature of strongly scat- 
tered (diffuson-like) waves in random media. 

The following picture Jl|-Q of vibrational modes in 
glasses has emerged from numerical studies of realistic 
models of vibrational disorder |Q-^,0-^o) . Figure |l| shows 
the calculated vibrational density of states for a model of 
a glass used in this paper (amorphous Sio.gGeo.i), with 
three specified spectral regions. At the lowest frequencies 
vibrations are acoustic-phonon-like propagons, weakly 
scattered sound waves with well defined momentum and 
polarization. At low frequencies one also finds resonance 
modes which have unusually large amplitude at certain 
inhomogeneous regions, while resembling propagons else- 
where. As the frequency increases, propagon momenta 
become less and less certain until the frequency reaches 
the so called Ioffe-Regel limit where the mean free path 
becomes comparable to the wavelength and the wave vec- 



tor concept is no longer valid. The emerging modes are 
diffusons-extended modes not able to propagate ballisti- 
cally, rather spreading out in a special diffusive fashion. 
Diffusons, forming the majority of the spectrum, are the 
most natural modes for glasses. Atoms in glasses display 
a short-range order by globally preserving local quantities 
like interatomic distances or coordination numbers. But 
at large distances the order is lost. Diffusons have the 
same property: displacements of neighboring atoms are 
strongly correlated, while those between distant atoms 
are uncorrelated and the global displacement pattern, 
like the structure of the glass, appears random. Finally, 
the highest-frequency modes are localized-locons. The 
purpose of this letter is to compare spectral properties of 
propagons and diffusons with random matrix theory and 
theories of wave- like modes in disordered systems, to sort 
out universal and nonuniversal features. 

Random matrix theory was originally devised to deal 
with the spectra of complicated nuclei, but by now it 
has been successfully applied in many areas of physics 
J5| . The main assumption is that certain spectral proper- 
ties of complex (as opposed to simple) Hamiltonians are 
universal, not dependent on the concrete realization of 
complexity; these properties can be therefore calculated, 
in many cases analytically, as averages over an ensem- 
ble of (similarly) complex Hamiltonians. I focus here on 
the Gaussian orthogonal ensembles (GOE) which are sets 
of real symmetric matrices with elements randomly se- 
lected from a Gaussian probability distribution. Dynam- 
ical matrices which determine the frequencies of atomic 
vibrations are also real symmetric so, if their structure is 
sufficiently complicated, their properties should be sim- 
ilar to the properties of the GOE. For the latter theory 
predicts that neighboring levels (frequencies) repel each 
other: there is a zero probability of finding two equal 
levels. A convenient measure of level repulsion is the 
level-spacing distribution (LSD) p which, for the GOE 
ensembles, is very accurately expressed by the Wigner 
surmise B 



Pw(s) = -sexp(--s ), 
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with s denoting the level spacing. Random matrix the- 
ory has also predictions about correlations between more 
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FIG. 1. Calculated vibrational density of states (shaded) 
for the 1000-atom model of amorphous alloy Sio.gGeo.i. 
The vertical lines are the Ioffe-Regel limit at 13 meV and 
the mobility edge at 70 meV, dividing the spectrum into 
propagons (below the Ioffe-Regel limit), diffusons (between 
the Ioffe-Regel limit and the mobility edge) , and locons (above 
the mobility edge). Also plotted is the inverse participation 
ratio 1/p (dashed) and the Thouless number Nt (solid) of the 
pure amorphous silicon model. 



distant levels. Such correlations can be quantified in sev- 
eral ways. Here I use the £2 statistics which measures, 
given an interval of size N, the variance of the number 
of levels in that interval. For GOE at large N 0] : 
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where 7 « 0.5772 is Euler's constant. The logarithmic 
dependence expresses the rigidity of a GOE spectrum. 

I evaluate LSD and £2 for the vibrational spectrum of 
a model of amorphous Sio.gGeo.i alloy whose atomic co- 
ordinates and interatomic forces are taken from the real- 
istic model of amorphous silicon @,[ll| generated by ran- 
domizing a crystalline silicon structure (Wooten-Winer- 
Weaire recipe [[l2|) and then relaxing the atoms to a lo- 
cal minimum of the Stillinger- Weber interatomic poten- 
tial [ fl3| |. The model structure and physical properties 
agree very well with experiment ||,[7|j8). Here I use a 
model with 1000 atoms arranged in a cube of side 27 A. 
The Sio.gGeo.i alloy is obtained by substituting germa- 
nium atomic masses for the masses of randomly picked 
10% silicon atoms. The mass disorder is small to qual- 
itatively change the vibrational properties of the origi- 
nal amorphous silicon model (similar amorphous silicon- 
germanium alloys were studied in Refs. |7 14|), but large 
enough to generate a useful ensamble. Below I use 
N r = 3500 different mass-disorder realizations of the al- 
loy. The model has 3000 vibrational frequencies whose 
spectrum is in Figjl] (averaged over the ensamble), along 
with the alloy inverse participation ratio 1/p (averaged 
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FIG. 2. Calculated level-spacing distribution (LSD) for 
propagons and diffusons, and the GOE result (shaded), Eq|j]. 

over 15 samples), allowing an identification of the mobil- 
ity edge at about 70 meV. The Ioffe-Regel limit is also 
indicated at 13 meV (taken from the studies of the pure 
silicon models |7j ; the actual transition is probably quite 
broad, perhaps plus and minus 2 meV). For the statis- 
tics, I take for propagons all the modes, about 200, from 
5 to 13 meV (this includes a small number, less than 5%, 
of resonance modes which behave mostly as propagons 
anayway) and for diffusons all the modes, about 1000, 
from 35 to 65 meV (which is enough away from the locons 
and from the anomalous quasi-localized region around 30 
meV at the DOS dip and 1/p peak). 

To calculate spectral averages that can be compared 
with the predictions of random matrix theory, the vi- 
brational spectrum has to be "unfolded," that is, locally 
rcscalcd to compensate for the spectral variations of the 
local average level spacing. A standard unfolding pro- 
cedure H involves finding the average cumulative DOS 
(a staircase function) of the spectral region of interest. 
The unfolded spectrum is then the set of values of the 
average function at the original frequencies. There are 
two ways of finding the average: by ensamble averag- 
ing or by analyzing a smooth fit to a single spectrum. 
The latter method works fine if an analytical behavior 
of the cumulative DOS is known. Without an analytical 
guidance, however, the procedure is more art than sci- 
ence JTq] . This is why I use amorphous alloy Sio.gGeo.i 
(rather than just pure amorphous silicon) for spectral 
statistics: first, there is no analytical formula available 
for DOS of amorphous silicon and, second, mass disor- 
der is much easier to implement than structural disorder. 
With a sufficiently large ensamble at hand, the averages 
are straightforward to evaluate. If p l c (oj) is the cumula- 
tive DOS of the z-th sample at frequency o>, the average 
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FIG. 3. Calculated E2 for propagons and diffusons. The 
dashed line is the GOE result and the thin line accompa- 
nying the propagon curve is E 2 = 0.0076 x(N/l.l)' s/2 , the 
Altshuler-Shklovskii power law, Eq. ^, with effective Thou- 
less number 1.1. The insets are the log- log plots for propagons 
(up) and diffusons (down). 

is p c (w) = (1/N r ) ^iP 1 ^)- After obtaining N r unfolded 
spectra, LSD and £2 for a certain frequency range are 
calculated by averaging over both the spectral region and 
the ensamble 

The calculation of LSD for both propagons and dif- 
fusons is in Fig. |^. The agreement with the Wigner 
surmise is excellent, proving the level repulsion between 
neighboring modes and, once again, the robustness of the 
LSD predictions. This connection between vibrations in 
glasses and GOE matrices was already reported in Ref. 
H (and in Ref. [|l7j for other disordered-lattice systems). 
It turns out, however, that random matrix theory is not 
obeyed by the correlations between more distant than 
neighboring (as in LSD) levels. 

Consider first propagons. The calculated £2 is in Fig. 
|§|. Deviations from random matrix theory, starting al- 
ready at small N, grow as a power law rather than a 
logarithm. The power law, which extends from N « 40 
to 200 is identified as £ 2 ~ N 3 / 2 (best seen in the in- 
set), and is precisely the one predicted by Altshuler and 
Shklovskii Q , originally for weakly scattered electrons on 
a 3-dimensional disordered lattice: 

£3 (AO ~ 0.0076 x (N/N T Y /2 ; (3) 

the exact proportionality constant is ^/2/Qtt 3 0.0076. 
The scaling parameter Nt is the Thouless number, Nt = 
HD/SL 2 , the uncertainty in the energy (in the units of 
the average local level spacing S) of a wave-packet-like 
state diffusing with diffusivity D throughout a sample of 
size L. Equation U is valid for N 3> Nt || and propagons 
obey it with N T « 1.1 (see Fig. Figure [l] plots N T 
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FIG. 4. Spectral statistics E2 of selected diffuson regions. 
The intervals in the legend and with the two of the curves 
are the spectral ranges in meV (see DOS in Fig.0). The GOE 
curve is also shown. The Thouless numbers of the regions are: 
3.2 (18-22), 1.2 (44-48), 1.3 (48-52), 1.2 (52-56), 0.96 (56-60), 
and 1 (60-64). The corresponding slopes of the linear increase 
are: 0.011, 0.049, 0.063, 0.085, 0.074, and 0.068. For example, 
diffusons from 44 to 48 meV have Nt ~ 1.2, and at large N 
their S 2 ~ 0.049 x N. 

as a function of vibrational frequency for the pure silicon 
1000-atom model, calculated from the model diffusivity 
D [Q. The result for the alloy would look identical, per- 
haps being smaller by a few percent (D for a similar alloy 
with 25% of substituted Ge atoms differs only slightly 
from D of the pure silicon sample M). We see that Nt 
is frequency dependent-it varies from below 0.1 to above 
4 in the propagon region. Considering that Eq. [| comes 
with some effective Nt, biased towards lower values, the 
result for the effective Nt ~ 1-1 is reasonable (the aver- 
age of 1 /N^! 2 over the propagon frequencies gives about 
1.7; for the alloy it will be slightly less), confirming the 
theory of Ref. §. 

Diffusons show nonuniversal behavior too. In Fig. |^, 
their £2 grows linearly with increasing TV for about a 
decade, from 30 to almost 400 (the growth slows down 
for greater N, most likely due to finite size effects). Like 
propagons, diffusons become less spectrally correlated 
when separated by more than a few other modes (of or- 
der 1; the deviation at small N is not visible because 
of the graph scale). But unlike propagons, the power 
law is linear, £2 ~ N. This different dependence is not 
surprising, as diffusons are not weakly scattered modes 
envisioned in the derivation of Eq. [| Diffusons do diffuse 
with time p8| , but in a different way than do propagons 
or electrons (in weakly disordered systems) . Equation ^ 
was derived by perturbation theory which assumes elec- 
tron states to be plane waves between scattering events. 
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In other words, the mean free path of electrons is much 
greater than their wavelength: the electrons are below 
the Ioffe-Regel limit and their diffusion resembles a ran- 
dom walk (like propagons). Diffusons, on the other hand, 
diffuse more like a free-particle quantum mechanical wave 
packet (this type of diffusion was branded "intrinsic" 0); 
diffusons are not perturbative states. 

To explore possible sources of differences between £ 2 
for propagons and diffusons, it is instructive to repeat 
the qualitative argument S which leads to Eq. 0. A 
given state can mix with neighboring N 3> Nt states 
if it diffuses a time tn ~ h/NS, or over a region with 
size L N = {Dtn) 1 ' 2 = (N T /N)^ 2 L. The N states in 
a volume (L^) 3 will be correlated as predicted by ran- 
dom matrix theory, with the uncertainty in the number 
of levels of order unity (a rigid spectrum). But states 
in different regions (all of size L/v) of the sample will 
be uncorrelated, so the uncertainty (measured by £ 2 ) in 
the number of states in the whole sample will be propor- 
tional to the number of the independent regions/spectra. 
This number is (L/L N ) 3 = (N/Nt) 1 ' 5 which is Eq. |. 
The argument can be made (naively) more general by 
considering an anomalous diffusive transport Ln ~ t^J 
(a = 1/2 in the usual case) and a fractal filling of volume 
lP (J3 = 3 for uniform filling). Then £ 2 ~ N a P, and 
for diffusons one gets a(3 = 1. If the space filling is or- 
dinary ((3 = 3), the diffusive transport is not (a = 1/3). 
If the transport is ordinary (a = 1/2), the spreading is 
fractal ((3 = 2). It is not clear wheather diffusons behave 
one way or the other (or both), or whether the above 
reasoning, which is a pure speculation, makes sense. 

If there is a scaling of the type £ 2 ~ (N/Nt), the best 
way to see it is to look at spectral statistics of small dif- 
fuson regions of almost constant Nt] the regions must 
contain enough levels to see the systematic deviations 
from GOE. Figure ^ plots E 2 of six such regions, all dis- 
playing the linear law S 2 ~ N, all with different propor- 
tionality coefficients (note that the law holds also for the 
18-22 region not considered in the calculation of Fig. ||). 
From the values of the coefficients and the regions' aver- 
age Nt's given in the caption, one cannot see any obvious 
relationship between the coefficients and Nt- There is 
certainly no support from the available data of a scaling 
of the type S 2 ~ (N/Nt)'. by plotting S 2 versus N/N T 
(not shown) the curves do not collapse onto or close to a 
single line. While there definitely is a dependence of E 2 
on Nt (explicit or implicit), the present spectrum size 
does not allow mc to decipher it and considering larger 
models is presently unfeasible. 

I thank P. B. Allen, J. Fcldman, and A. M. Halasz for 
stimulating discussions. 
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